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Abstract 

The present results on the — » v T oscillation search by the CHORUS experiment at CERN are summarised. A 
fraction of the neutrino interactions collected in 1994-1995-1996 by the CHORUS experiment has been analysed, 
searching for v T charged current interactions followed by the r lepton decay into a negative hadron or into a muon. 
A sample of 68, 156 events with an identified final state muon and 7,206 events without an identified muon in the 
final state have been located in the emulsion target. Within the applied cuts, no v T candidate has been found. This 
result leads to a 90% C.L. limit — » u T ) < 6.0 ■ 10~ 4 on the mixing probability. 
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The CHORUS experiment has recently reported [jj], 
^[ [jj a limit on — v T oscillation obtained by the anal- 
ysis of a subsample of neutrino interactions, taken in 
1994-1995, both with and without an identified pT in 
the final state. This paper contains an update of the statis- 
tics of the subsample of events with an identified pT in 
the final state. 

The experimental setup and the characteristics of 
the CERN wide band neutrino beam are summarised 
in |[TJ and described in more details in [|J. 

1 The apparatus 

The " hybrid" CHORUS apparatus combines a 770 
kg nuclear emulsion target with various electronic de- 
tectors : a scintillating fibre tracker system, trigger ho- 
doscopes, an air-core magnet, a lead/scintillator calorime- 
ter and a muon spectrometer. 

Details about the experimental setup and the perfor- 
mances of the sub-detectors can be found in [Q[. 

2 Data collection and event selection 

In the 1994-1997 period, CHORUS has collected 
2,271,000 triggers corresponding to 5.06 10 19 protons 
on target. Of these, 458,601 have a muon identified in 
the final state (the so called l/i events) and 116,049 do 
not (the so called 0/i events) and a vertex position com- 
patible with one of the four emulsion target stacks. 

All tracks, associated to the interaction vertex, with 
an angle less than 0.4 rad with the beam axis and, in 
view of the large emulsion background of muons origi- 
nating from a nearby secondary target, bigger than 0.05 
rad from the direction of this target are extrapolated down- 
stream and selected for further analysis. These tracks are 
searched for in the emulsion if their charge is negative 
and their momentum is in the range < ph < 20 Ge V/c 
and < p p < 30 Ge V/c for hadrons and muons, respec- 
tively. 

It should be noted that use has not been made of the 
energy deposition in the calorimeter to reject electrons 
or unidentified muons. The contribution of the r~ lep- 
tonic decay modes to the 0/i data sample is taken into 
account in the evaluation of the sensitivity. 

3 Scanning procedure 
3.1 Vertex location 

The various steps leading to the plate containing 
the vertex by means of fully automatised microscopes 
are identical to those described in [jj], pp. They are inde- 
pendently applied to all the selected tracks in the event, 
the muon for the l/i events, all the negative tracks for 
the 0/i events. A track which has been found in the in- 
terface emulsion sheets is followed upstream in the tar- 
get emulsion, using track segments reconstructed in the 
most upstream 100 /im of each plate, until the track dis- 
appears. This plate is referred to as the vertex plate, 
since it should contain the primary neutrino vertex or 
the secondary (decay) vertex from which the track orig- 
inates. The three most downstream plates of each stack 
are used to validate the matching with the interface emul- 
sion sheets and are not considered as possible vertex 
plates. The scanning results are summarised in Table [j]. 



The mean efficiency of this scan-back procedure is 
found to be ~ 32% and ~ 42% for 0/i and l/i events, 
respectively. A detailed simulation of the scanning cri- 
teria shows that the difference mainly reflects the poorer 
quality of the hadron track predictions at the interface 
emulsion sheets, because of the difficulty of reconstruc- 
tion inside a dense hadronic shower and the larger mul- 
tiple scattering owing to the lower average momentum. 

3.2 Decay search 

Once the vertex plate is defined, automatic micro- 
scope measurements are performed to select the events 
potentially containing a decay topology (kink). Differ- 
ent algorithms have been applied, as a result of the pro- 
gress in the scanning procedures and of the improving 
performance in speed of the scanning devices. They are 
described in [[[], pj] and briefly recalled here. 

In the first procedure the event is selected either 
when the scan-back track has a significant impact pa- 
rameter with respect to the other predicted tracks or when 
the change in the scan-back track direction between the 
vertex plate and the exit from the emulsion corresponds 
to an apparent transverse momentum, px, larger than 
250 Me V/c. For the selected events and for those with 
only one predicted track, digital images of the vertex 
plate are recorded and are analysed off-line for the pres- 
ence of a kink. 

The second procedure is restricted to the search of 
long decay paths. In that case the vertex plate is assumed 
to contain the decay vertex of a charged parent particle 
produced in a more upstream plate. With this procedure 
only kink angles larger than 0.025 rad are detected. 

For the events selected by either one of these pro- 
cedures, a computer assisted eye-scan is performed to 
assess the presence of a secondary vertex and measure 
accurately its topology. A decay candidate must sat- 
isfy the following criteria: 

1. the secondary vertex appears as a kink without 
black prongs, nuclear recoils, blobs or Auger elec- 
trons; 

2. the transverse momentum of the decay muon (ha- 
dron) with respect to the parent direction is larger 
than 250 Me V/c (to eliminate decays of strange 
particles); 

3. the kink, in the 0/i channel, occurs within 3 plates 
downstream from the neutrino interaction vertex 
plate. Because of the lower background, the kink 
search in the muonic decay channel was extended 
to 5 plates, with a gain in efficiency of about 8%. 

No r~ decay candidate has been found satisfying 
the selection criteria. 

4 Experimental check of the kink finding 
efficiency 

The kink finding efficiency, ekink, has been eva- 
luated by Monte Carlo simulation and experimentally 
checked looking at hadron interactions and dimuon events. 

A sample of about 55 m of hadron tracks was scanned 
in emulsion. In the decay search procedure, 21 neutrino 
interaction events with a hadron interaction have been 



Table 1: Current status of the CHORUS analysis 







1994 


1995 


1996 


Protons on target 




0.81 • 10 19 


1.20 • 10 19 


1.38- 10 19 


Emulsion triggers 




422,000 


547,000 


617,000 


l/i to be scanned 




66,911 


110,916 


129,669 


O/i to be scanned 




17,731 


27,841 


32,548 


l/i scanned so far 




42,154 


49,912 


72,615 


O/i scanned so far 




8,908 


12,635 




l/i vertex located (in 


the 33 most upstream plates) 


18,286 


20,642 


30,128 


O/i vertex located (in 


the 33 most upstream plates) 


3,401 


3,805 






detected. This result is in good agreement with the ex- 
pected value of (24 ± 2) from Monte Carlo simulation. 

Part of the dimuon sample (two muons on the fi- 
nal state) collected in 1995 and 1996 have been anal- 
ysed searching for the reaction v^N — > /i~ D + X with 
the subsequent muonic decay of the D + . Assuming a 
charm yield of about u c harm/ ccc ~ 5%, in the sam- 
ple scanned we expect (22.8 ± 3.9) dimuon events and 
we found 25 events. This result shows that Monte Carlo 
simulation and real data are in good agreement. 

Although the number of events is too low to draw 
quantitative conclusions, we can take these results as 
qualitative checks of the simulation of the automatic scan- 
ning procedure. 

5 Sensitivity and backgrounds 

In this section we discuss the expected background 
from known sources and, in absence of t~ candidates, 
limits on the — > v T oscillation parameters are de- 
rived. Both signal efficiencies and background estima- 
tion have been evaluated from large samples of events, 
generated according to the relevant processes, passed 
through a GEANT [||] based simulation of the detec- 
tor response. The output is then processed through the 
same reconstruction chain as used for data. The simu- 
lated tracks in emulsion are used to estimate the effi- 
ciencies of each scanning step. Apart from the kink de- 
tection efficiency, only ratios of acceptances enter the 
calculation of the experimental result and most of the 
systematic uncertainties of the simulation cancel out. 

5.1 Background estimates 

Sources of background for the hadronic r decay 
channel are: 

- the production of negative charmed particles from 
the anti-neutrino components of the beam. These 
events constitute a background if the primary //+ 
or e + remains unidentified. Taking into account 
the appropriate cross-sections and the branching 
ratios, in the present sample we expect ~ 0.02 
events from these sources; 

- the production of positive charmed mesons in char- 
ged current interactions, if the primary lepton is 



not identified and the charge of the charmed par- 
ticle daughter is incorrectly measured. We expect 
~ 0.03 events from this source in the present sam- 
ple; 

- the associated charm production both in charged 
(when the primary muon is lost) and neutral cur- 
rent interactions, when one of the charmed parti- 
cles is not detected. The cross-section for charm- 
anticharm pair production in neutral current inter- 
actions relative to the total charged current cross- 
section has been measured by the E531 experi- 
ment [§] to be 0.13l° 11%. The production rate 
of associated charm in charged current interac- 
tions is unknown, but an upper limit is available 
(< 0.12% In the present sample, the esti- 
mated background from this process represents 
< 0.01 events. 

- the main potential background to the hadronic t~ 
decays is due to so-called hadronic "white kin- 
ks", defined as 1 -prong nuclear interactions with 
no heavily ionising tracks (black and grey tracks 
in emulsion terminology) and no evidence for nu- 
clear break up (evaporation tracks, recoils, blobs 
or Auger electrons). Published data allowing to 
determine the white kink interaction cross-section 
are scarce [|7| |J]. In a dedicated experiment with 
4 GeV pions at KEK [Q], a very steep fall-off in 
Pt was observed and only 1 out of 58 observed 
kinks had a pr larger than 300 MeV/c. 

Since the experimental information of the pt de- 
pendence is statistically poor at large values, a 
Monte Carlo simulation, based on a modified ver- 
sion of FLUKA [^, [lot ], has been performed. The 
results of this simulation are in good agreement 
with the pt dependence of the KEK measurement. 
An effective white kink mean free path in emul- 
sion, X w k ~ 22 m, is obtained for a pr cut at 250 
MeV/c, using a pion energy spectrum as observed 
in the 0-/i sample. The above result is compatible 
with the observation of 4 events with pr > 250 
MeV/c, all at large distance from the primary ver- 
tex, along a total path of ~ 92 m of scan-back 
tracks, and corresponds to a background estimate 



2 



of 0.5 events within 3 plates downstream from 
the primary vertex plate. As the statistics will in- 
crease, this background source will be better mea- 
sured using the CHORUS data. It will be possible 
to significantly reduce, by kinematical analysis, 
the effect of this background on the sensitivity to 
a possible oscillation signal. 
The main source of potential background in the muonic 
t channel is the charm production. We expect less than 
0. 1 events in the current sample from the anti-neutrino 
components of the beam: 



BR 



(1 or 4] 



= BR(t — > v T u,ji 



(17.35 ± 



t (D e )M -^ f i+(e+)D-X 



0.10)% [|2p). 

- BR 2 = BR(t -> v T h-nh°) = (49.78±0.17)%jPl); 

- BR 3 = BR(t -> v T D e e~) = (17.83±0.08)% 

- $> v the incident flux spectrum; 

- <7 T the charged current v T interaction cross-section; 

- A T i the acceptance and reconstruction efficiency 
for the considered channel (up to the vertex plate 
location); 

- t T i the corresponding efficiency of the decay search 
procedure; 

With proper averaging (denoted by ()), N Ti can also 
be written as a function of n,-: 



followed by 



D~ 



N Tl = BR, ■ n % • (P MT ) 



(gr) {Ki) 



(2) 



in which the fi + (e + ) escapes the detection or is not 
identified. 

The prompt v T contamination of the beam jTT[ | is 
a background common to both the hadronic and muonic 
decay channels. For the present sample the expected back- 
ground is much less than 0.1 events. 

5.2 Oscillation sensitivity 

In the usual approximation of a two-flavour mixing 
scheme, the probability of v T appearance in an initially 
pure beam can be expressed as 



where 



m = N, (the number of located charged current 
interactions corresponding to the considered 
event sample) and n% — na = = (N,) Q _ fi 
(the product of N, and the relative fraction of the 
0-fjL sample for which the analysis has been com- 
pleted); 

-dE. It takes into account 



da 



m(t) 



(°m(t)> = J 

quasi-elastic interactions, resonance production and 
deep inelastic interactions {<j{^-^) sys t ~ 7%); 



dE 



3>„ -dE 



1.27 



■sin 



V(E,L) ■ 
AmUeV 2 ) 



L(km) 



E(GeV) 



dL 



where 

- E is the incident neutrino energy; 

- L is the neutrino flight length to the detector; 

- 9, T is the effective — v T mixing angle; 

- Am^ T is the difference of the squared masses of 
the two assumed mass eigenstates; 

- ty(E, L) is the fraction of with energy E orig- 
inating at a distance between L and L + dL from 
the emulsion target. 

The t~ channels considered in the — > v r oscil- 
lation search we describe in this paper are: 

1) t — » /j,, 2) t — > h, 3) r — > e and 4) t — ► p, (the /i 
is not identified) channels. 

The expected number, N T i (i = 1, 2, 3, 4), of ob- 
served r~ decays into a channel of branching ratio BRi 
is then given by 



dE 

syst ~ 7%); 

(e T ,) is the average efficiency of the decay search 
procedure for the accepted events (cr((e T i)) sys t ~ 
10%); 

To allow an easy combination of the results from 
the l-fi and 0-u event samples, it is useful to define the 
"equivalent number of muonic events" of the 0-/x sample 
by 



E 



{A Ti ) {e Tl ) BR, 



(A Tfi ) (e rM ) BR, 



(3) 



The 90% C.L. upper limit on the oscillation proba- 
bility then simplifies to 



P 



< 



2.38 -twa 



where r„ 



BR, ■ (e T ,) ■ [N M + N^} 
(cr M )/(cr T } and r A = (A,)/ {A TfJ _). 



(4) 



with 



BRi 



• Put ■ o> ■ A Tl ■ e Ti ■ dE (1) 



In the above formula, the numerical factor 2.38 takes 
into account the total systematic error (17%) following 
the prescription given in jl3|]. The systematic error is 
mainly due to the reliability of the Monte Carlo simula- 
tion of the scanning procedures. 

The estimated values of the quantities appearing in 
this expression are given in Table |[ No statistical errors 
are quoted since they are much smaller than the system- 
atic uncertainty. 
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12] 
[3] 



14] 
[5] 
16] 
17] 
18] 
19] 



sin 2 ?^ 

Figure 1 : Present result compared with the recent NO- 
MAD result [ 14 1 (full line) and the previous limits (dot- 
ted lines). 



Using the present sample the following 90% C.L. [10] 
limit is obtained 



P^t < 6.0 ■ 10" 



(5) 



In a two flavour mixing scheme, the 90% C.L. ex- 
cluded region in the (sin 2 26^ T , Am^ T ) parameter space 
is shown in Figure [j]. Maximum mixing between and \\\\ 



v T is excluded at 90% C.L. for A m 2 T > 0.9 eV 2 ; the 
large Am 2 are excluded at 90% C.L. for sin 2 29 tlT > 
1.2 • 10~ 3 . 



[12] 
[13] 

[14] 



6 Conclusions 

The emulsion scanning methods, previously described 
in [|lj |||], have been applied to a fraction of the 1994- 
1995-1996 data. No r~ decay candidate has been found, 
leading to a more stringent 90% C.L. upper limit on the 
Vjj, — > v T oscillation probability (P^ T < 6.0-10 -4 ). This [15] 
negative result is compatible with an estimated back- 
ground of less than one event, mainly from "white kink" 
secondary interactions in the 0/i channel. With the large 
increase in statistics expected from the ongoing analy- 
sis of the data, a direct measurement of this background 
process will be possible and kinematical cuts are planned 
for its reduction. Furthermore, a second phase of the 
analysis (with better efficiencies, larger statistics and faster 
automatic emulsion scanning) has started with the aim 



of reaching the design sensitivity (P MT < 1.0-10 )[15| 
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Table 2: Quantities used in the estimation of the sensitivity 





1994 


1995 


1996 




18,286 


20,642 


30,128 


r<? 


1.89 


1.89 


1.89 




0.93 


0.93 


0.93 




0.39 


0.39 


0.39 


(A Th ) 


0.17 


0.17 




{A Te ) 


0.093 


0.093 




{Am) 


0.026 


0.026 




{trii) 


0.53 


0.35 


0.37 


(erh) 


0.24 


0.25 




(ere) 


0.12 


0.13 




{e T p) 


0.22 


0.23 




N eq 

u 


11,987 


12,743 





